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Abstract 

Deep inelastic ep scattering data, taken with the HI detector at HERA, are used to 
study the event shape variables thrust, jet broadening and jet mass in the current 
hemisphere of the Breit frame over a large range of momentum transfers Q between 
7 GeV and 100 GeV. The data are compared with results from e+e“ experiments. 
Using second order QGD calculations and an approach to relate hadronisation effects 
to power corrections an analysis of the Q dependences of the means of the event 
shape parameters is presented, from which both the power corrections and the 
strong coupling constant are determined without any assumption on fragmentation 
models. The power corrections of all event shape variables investigated follow a 
1/Q behaviour and can be described by a common parameter do. 
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1 Introduction 


Hadronic final states in deep inelastic ep —>■ e'X scattering (DIS neutral current interaction) 
offer an ideal environment to study hadronisation phenomena and to measure the strong coup¬ 
ling constant over a wide range of momentum transfer Q in a single experiment. Event shape 
variables have been investigated in experiments and used to extract the strong coupling 
constant as{Mz) independent of any jet algorithm, see e.g. ref. [^j. In deep inelastic scattering 
a similar analysis can be performed, provided the current (quark) fragmentation can be iso¬ 
lated from the target (proton remnant) fragmentation region. A particularly suitable frame of 
reference in which to study the current region with minimal contamination from target frag¬ 
mentation effects is the Breit frame. Consider, for illustration, the quark parton model. In the 
Breit system the purely space-like gauge boson 7 /Z with four-momentum q^/z = {0, 0 , 0 , —Q} 
collides with the incoming quark with longitudinal momentum p*” = Q/2. The outgoing quark 
is back-scattered into the current hemisphere with longitudinal momentum p°“* = —Q/2, while 
the proton fragments in the opposite hemisphere. The available energy for fragmentation in the 
Breit current hemisphere is Q/2. This situation can be compared with qq annihilation, 

where the available energy in one hemisphere, e.g. dehned by the thrust axis, is ^/s^/2. DIS 
event properties in the Breit current hemisphere should be similar and may be comparable to 
those in e’''e“ collisions at a scale Q = ^ys^. 

In this paper an analysis of the event shape variables thrust, jet broadening and jet mass 
in the current hemisphere of the Breit frame in deep inelastic ep scattering is presented 
and compared with results from e+e“ experiments. The kinematic phase space extends over 
7 GeV < Q < 100 GeV in momentum transfer and 0.05 < y < 0.8 in the kinematic variable y. 
Based on only recently available second order QGD calculations 0.0 and on new developments 
in the treatment of non-perturbative hadronisation contributions to hnal states in deep inelastic 
scattering the Q dependences of the mean event shape parameters are investigated to study 
simultaneously both the power, or hadronisation, corrections and the strong coupling constant 
as{Mz) independently of any fragmentation model. 


2 Data Selection and Analysis Procedure 

The HI Detector Deep inelastic ep scattering events were collected with the HI detector 
at HERA {Ef. = 27.5 GeV, Ep = 820 GeV, ^/s = 300 GeV). The HI detector is described 
elsewhere [Q; here only the components relevant for the present analysis are briefly recalled. 
Only calorimetric information is used to reconstruct electromagnetic and hadronic clusters. The 
event vertex and the direction of the scattered lepton are obtained with the tracking detectors. 

The calorimeter system consists of a liquid argon (LAr) calorimeter, a backward calorimeter 
and a tail catcher (instrumented iron yoke). The LAr sampling calorimeter covers the polar 
angleQ range 4° < 0 < 153° and all azimuthal angles. It consists of a lead/argon electromagnetic 
section followed by a stainless steel/argon section for the measurement of hadronic energy. 
Electromagnetic energies are measured with a resolution of a{E)/E ~ 12%/\fE © 1% and 
hadronic energies with a{E)/E ~ 50%/\/E © 2%, as obtained from test beam results. The 
absolute energy scales are known to 3% and 4% for electrons and hadrons, respectively. A 
lead/scintillator electromagnetic backward calorimeter (BEMG) extends the coverage at large 

^Polar angles 9 are defined with respect to the incident proton direction. 
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angles (155° < 0 < 176°). The electromagnetic energy calibration is known to 1%. Since 
1995 the backward region has been eqnipped with a lead/scintillating hbre calorimeter The 
instrnmented iron flnx retnrn yoke is used to measure the leakage of hadronic showers. 

Located inside the calorimeters is a tracking system, which consists of central drift and 
proportional chambers (25° < 6 < 155°), a forward track detector {7° < 6 < 25°) and a 
backward proportional chamber (155° <d< 175°). In 1995 the latter was replaced by backward 
drift chambers. The tracking chambers and calorimeters are surrounded by a superconducting 
solenoid providing a uniform held of 1.15 T throughout the tracking volume. 

The luminosity is determined from the rate of the Bethe-Heitler process ep —> epy mea¬ 
sured in a luminosity monitor far downstream in the electron direction. 


Data Selection The DIS data cover a large range of momentum transfer Q between 7 GeV 
and 100 GeV and are selected with the following criteria: 

1. The energy of the isolated scattered lepton has to exceed 77' > 10 GeV. 

2 . The polar angle of the scattered lepton has to be within 157° < 9e < 173° (low Q sample) 
or 30° < Of. < 150° (high Q sample). 

3. The hadronic energy clusters have to be well contained within the calorimeter acceptances 
5.7° <eh< 170°. 

4. There must be hadronic energy in the forward region (4° < 6* < 15°) which is larger than 
0.5 GeV in order to exclude diffractive events with large rapidity gaps. 

5. The energy and longitudinal momenta must satisfy 30 GeV < Ei (1—cos Oi) < 65 GeV, 

where the sum extends over all energy clusters. 

6 . The total hadronic energy in the Breit current hemisphere has to exceed 0.1 Q. 

7. The kinematic variables y have to fulhll 0.05 < ye < 0.80 (measured using the lepton) 
and 0.05 < yh (measured using the hadronic energy flow). 

For the low Q event sample [Q = 7—10 GeV, 1994 data with beam, integrated luminosity 
C ~ 2.9 pb~^) the lepton is detected in the BEMG. The lepton direction is measured using the 
event vertex and the backward proportional chamber in front of the calorimeter. For the high 
Q event sample (Q = 14 — 100 GeV, 1994 - 1996 data with beams, C ~ 11.7 pb“^) the 
lepton is detected in the LAr calorimeter and its direction is taken from an associated track 
measured in the central tracking system. 

The event selection criteria are chosen so as to ensure a good measurement of the hnal state 
and to provide a clean DIS data sample. Requirements 1, 2 and 7 assure that the kinematic 
quantities are well measured. The criterion 5 provides a good containment of the hnal state 
particles and rejects events with a hard photon radiated from the incoming lepton and together 
with the high y cut of requirement 7 suppresses photoproduction events with a misidentihed 
lepton. Using tagged events, the photoproduction background is estimated to be about 3% in 
the low Q sample and to be neglible at high values of Q. To dehne shape variables, events are 
required to have hadronic activity in the Breit current hemisphere according to criterion 6. 
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Analysis Procedure The kinematic quantities needed to perform the Breit frame transfor¬ 
mation, the negative squared momentum transfer and Bjorken x = Q^jys^ are calculated 
from the scattered lepton (T^', 6*e, 0e) and hadron 0^, <ph) measurements using 



= JEeE^cos^l, 



(1) 


E' , Oe 



(2) 

y = Ve 

= 1 —sm -- 
Ee 2 

for 

ye > 0.15 , 

III 

s- 

_ (1 - COS^/j) 

2Ee 

for 

ye < 0.15 . 

(3) 


The momentum transfer is always best measured using the lepton. The kinematic variable 
y is taken from the lepton for sufficiently large values. However, since the resolution in r/e 
degrades severely at low values, y is determined using y^ if ye < 0.15. This procedure ensures 
least uncertainty in the Lorentz transformation to the Breit frame. 

The data are corrected for detector effects (acceptance, resolution, energy scale) using a full 
Monte Carlo simulation of the detector response. The Monte Carlo event selection and analysis 
are the same as for the data. The LEPTO event generator including matrix elements plus 
parton showers is used over the full Q range. At low Q, however, events generated with ARI¬ 
ADNE according to the colour dipole model are also included in order to enlarge the Monte 
Carlo statistics. Both simulations are based on the MRS H parton density parametrisations 
and they describe the data well. 

The observed differential event shape distributions {T)} are corrected with a matrix method 
pcor _ Cij In principle the matrix C may be obtained by an inversion of the matrix 

M. which transforms the ‘true’ values into observed values Mij including the 

detector effects as described by a Monte Carlo simulation. In practice, however, the inversion 
of M. poses mathematical problems and leads to instabilities and oscillations, unless extremely 
large Monte Carlo statistics are available. Instead, the probability densities correlating 
the observed and ‘true’ quantities, i.e. are used to get the correction 

matrix Cij = pfj^ / Y.k ■ This method may still have some model dependence, which could 
be overcome by an iterative procedure. However, as the data are described acceptably by the 
Monte Carlo simulations no iteration is performed. 


Systematic uncertainties in the data analysis are evaluated for the mean values of the event 
shape variables, which are computed from the differential distributions. The following sources 
are considered: (i) the influence of the unfolding procedure, i.e. applying the matrix method 
or a bin-by-bin correction method; (ii) a variation of the energy scales, which affects the Breit 
frame transformation but has a small effect on the event shape parameters, since they involve 
ratios of hadronic momenta or energies; and (iii) the use of the HERWIG event generator IJT^ 


with fragmentation properties which are different from LEPTO. All contributions are of similar 
size over the whole range of Q. QED radiative corrections, which are studied with the DJANGO 
event generator |H|, are found to have a negligible effect on the event shapes. The sum of all 
systematic errors of the mean event shape variables varies between ~ 2% at low Q values and 
up to ~ 8% at high Q values. 
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3 Event Shapes in the Breit Current Hemisphere 


The infrared safe event shape variables thrust Tc and T^, the jet broadening Be and the jet mass 
Pc are studied in the Breit current hemisphere. The dehnitions of the event shape variables 
are given below, where the sums extend over all hadrons h (being a calorimetric cluster in the 
detector or a parton in the QCD calculations) with four-momentum ph = p/i} fulhlling 
cos (p/i ■ n) >0. The current hemisphere axis n = {0, 0, —1} coincides with the virtual boson 
direction. 


Thrust Tn 


Tr = max 


Eh\Ph- ht 

E/i I Ph I 


Thrust T, 


T, = 


Eh IPh■ n I Eh I P 


2 h 


Eh I Ph I 

Jet Broadening Be 

Eh I Ph X n 


Be = 


2 Eh I Ph 


EhIPh 


Eh I P±h I 

2 Eh I Ph 


Jet Mass pc 


pc = 


m 2 

'W 


EhPh 


ut = thrust axis , (4) 


n = hemisphere axis , (5) 


n = hemisphere axis , (6) 


(7) 


The normalized differential spectra of the event shape variables in the current hemisphere 
of the Breit frame and their mean values as a function of Q are presented in hgs. |l| and |^. The 
mean values of the event shape parameters are also listed in table |^. They extend over a large 
range of momentum transfers Q between 7 GeV and 100 GeV. The corrected data, spectra and 
mean values, are well described by the LEPTO Monte Garlo generator for all Q. 

A common characteristic of the event shape spectra of 1 — Tc, 1 — T^, Be and pc is that 
they get narrower and more peaked towards lower values as Q, or the available energy in the 
Breit current hemisphere, increases. This means that the energy flow becomes more collimated 
along the event shape axis. This fact is also evident from the mean values of the event shape 
variables, which exhibit a clear decrease with rising Q. 

The energy dependence of the mean thrust and jet masses of the HI DIS ep data may be 
compared with available results 

event shape parameters in e p scattering are calculated in a single hemisphere, i. e. the current 
hemisphere of the Breit system, and that a ‘natural’ axis comparable with the virtual boson 
direction does not exist in e+e“ annihilation. There, the virtual boson is at rest and the event 
axis is typically chosen to be the thrust axis. 

The Q dependence of the means (1 — in Eg. ^ is in qualitative agreement with the 
e’''e“ measurements. However, the e^e“ data points, coming from many experiments, tend 


12 1 at the scale Q = Recall, however, that the 
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to lie systematically above the ep data at low values of Q < 15 GeV. This can be partially 
understood from the fact that in e+e” experiments thrust is calculated for the whole event and 
not in a single hemisphere, which tends to deplete the high thrust region. This property is also 
demonstrated by comparing the thrust spectra of £g. |l|a with the corresponding distributions at 
similar energies of the PLUTO and TASSO experiments |^. Furthermore, from the dehnition 
of eq. (^), it can be seen that thrust measures a mixture of the longitudinal and transverse 
momentum components with respect to the boson direction. In strongly collimated events the 
thrust axis is close to the 7 /Z axis and is essentially given by the longitudinal momenta. 
However, for more isotropic events, which frequently occur at low Q, the thrust axis flips by up 
to 90° and is dominated by transverse momenta, resulting in higher thrust values and thus 
mimicking jet-like event configurations. Such a situation does not occur in analyses. The 
thrust dehnition of eq. (||) is expected to be closer to the case. Indeed, with a proper 
rescaling in order to cover the same thrust range, the mean values (1 — Tz)/2 exhibit a stronger 
Q dependence, see £g. |^, and are comparable to (1 — Tee) at Q < 50 GeV. 


The mean jet masses (pc) in hg- 0d show a weaker Q dependence for ep scattering than for 
e’''e“ collisions, where the average of the heavy and light jet masses of both event hemispheres 
are plotted. The ep data are systematically higher at larger values of Q > 20 GeV. A possible 
explanation is the normalisation of Pc to in eq. (|^). In the Breit current hemisphere multi- 
parton hnal states tend to produce visible energy > Q/2, which is different from the quark 
parton model expectation and from the situation in e'''e“ annihilation. Indeed, a normalisation 
to instead of gives better agreement between ep and data, particularly at high 

values of Q. 


As argued above, one does not expect the event shapes distribution in deep inelastic ep 
scattering to correspond exactly to those for e’''e“ annihilation. There are other differences, 
which are not related to the analysis method. The e’''e“ data contain a considerably larger 
fraction of bottom quarks with different fragmentation properties than the light quarks, which 
may be important close to threshold energies. Furthermore, there are additional 0(as) processes 
in deep inelastic scattering. Besides the common final state gluon radiation off quarks, in DIS 
there are contributions from initial state gluon radiation and boson gluon fusion. 


4 QCD Calculations and Power Corrections 


Theoretical Framework Any Q or energy dependence of the event shape variables can have 
the following origins: (i) the logarithmic change of the strong coupling constant as{Q) oc 1/ In Q, 
and (ii) power corrections or hadronisation effects, which are expected to behave like 1/Q. 
Recent theoretical developments suggest that 1/Q corrections are not necessarily related to 
hadronisation, but may instead be a universal soft gluon phenomenon associated with the 
behaviour of the running coupling at small momentum scales [^. ‘Universal’ means that 


they could be expressible in terms of a few non-perturbative parameters with calculable process- 
dependent coefficients. The event shape data will now be analyzed by applying this approach 
which relates hadronisation effects in the hnal state observables to power corrections. 

The mean value of any infrared safe event shape variable {F) at a scale pn, taken to be 
fJ'R = Q, can be written according to ref. ^ as 


(F) = (F)P"''* + (F)P°" 


( 8 ) 






^^^pert ^ asi^R) + (*^2 + ^ In ^ Ci^ , 


(F)P°"' = ap^—ln^ 


16 Hi 

3 71 Hr 


fJ'R 

Hi 


aoini) - asinn) - ^ + ^ + 

ZTT \ Hi /Jq I 


( 9 ) 

. ( 10 ) 


where I3q = 11 — 2/3 Nf, K = 67/6 — 7 r ^/2 — 5/9 Nf and Nf = 5 flavours. 

The perturbative part, eq. (||), represents the second order QCD prediction, where the 
coefficients ci and C 2 are calculated in the MS scheme at the scale hr = Q- The power 


corrections, eq. (pTf ) with coefficients ap and p (p = 0 , except p = 1 for (Be)) depending 
on the observable F, contain a free non-perturbative parameter q.q{hi) fo be evaluated at 
some ‘infrared matching’ scale Aqcd ^ Hi ^ ARi conventionally chosen to be p/ = 2 GeV. 
The parameter ao{Hi) can be interpreted as an effective strong coupling constant below the 
matching scale hi, calculated in leading order. The dependence on hi is partially compensated 
by the hi dependence of the other terms in eq. (0)- The dependence on the renormalisation 
scale Hr should help to compensate the scale dependence of the perturbative part. The power 
corrections are expected to behave like 1/hr (he. p = 0 ) for the event shape variables discussed 
in this paper, except for the jet broadening which is supposed to have an additional ln(/i^/ hi) 
term (he. p = 1). However, the prediction for (B/) is considered to be less reliable Q. In 
this approach the data can be directly used to determine in a simultaneous fit the power 
correction parameters do(h/) and the strong coupling constant as{Mz) without assuming any 
fragmentation model. 


Perturbative Predictions In deep inelastic scattering the perturbative part, eq. (P), may 
be obtained via 

pen ^ lo-F^dF ^ ^ 

da 

To get (F)P®’’* to 0{al) the Taylor expansion of eq. (|^) shows that the integral in the numerator 
has to be evaluated in second order QCD, whereas the total cross section citot needs only be 
known to first order, because the numerator vanishes, F = 0 , in the quark parton model. 

The perturbative calculations can be performed with the MEPJET [0 and DISENT 
programs. Both programs treat deep inelastic ep scattering to 0{al) via single 7 exchange. 
The neglect of Z exchange has no influence on event shape parameters in the range of Q 
investigated. The lack of virtual two-loop diagrams, however, leads to the restriction that extot 
can only be calculated to 0{as), which, as explained above, is sufficient for the present analysis. 

MEPJET applies the so-called ‘phase space slicing method’ for the integration. The real 
emission of partons (quarks, antiquarks and gluons) is calculated exactly for a two-parton 
resolution parameter Sij > Smin, the invariant mass squared, with typical values of Smin ~ 
0.01 GeV^. Soft and collinear approximations are used in the region where at least one parton 
pair has Sij < Smin- The final state infrared and collinear divergences cancel against virtual loop 
and tree diagrams. The cut on prevents the integration of the event shape distributions 
F da/dF being carried out over the whole phase space. Instead, a lower bound F^nt > 0, which 
depends on the event shape under consideration, has to be imposed. 

DISENT uses the ‘subtraction method’ and applies dipole factorisation formulae for the 
same set of 0(ag) diagrams. Since no Smin cut is needed and special care has been taken 
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concerning the numerical integration of F da/dF, which still contains integrable singularities, 
one is allowed to use the complete phase space 0 < F < F^ax- Therefore the DISENT results are 
applied throughout the present analysis. In the overlapping phase space regions with F > F^ut 
the differential event shape distributions of MEPJET and DISENT agree very well and the 
mean values from both programs yield compatible results to within the statistical accuracy, 
which is less than 1% in Erst order and about 3% in second order perturbative QCD. 


Fig. I shows the differential event shape spectra at the parton level based on the second 
order QCD calculations in comparison with the data. The lepton and partons are selected in 
the same way as in the data by requiring the criteria 1, 2, 6 and 7 of section Q. One clearly 
sees the strong impact of hadronisation effects, particularly at low values of Q. Towards higher 
values of Q the hadronisation becomes less important and the data and parton distributions 
approach each other. It is a consequence of the concept of ref. [Q that the mean values of the 
event shape spectra of the observed data and of the partons can be related by a single power 
correction parameter do- 


For each event shape variable the coefficients ci and C 2 are obtained from a £t of eq. to the 
Q dependence of the mean values, which themselves are calculated with DISENT from eq. (0. 
The theoretical calculations are performed by using the MRS H parton distributions 0 and the 
corresponding value of the strong coupling constant. The coefficients ci, C 2 of the perturbative 
part, eq. (H), and the coefficient @] of the power corrections, eq. (pRl), which are used in the 
following QCD analysis of the event shape variables are given in table 0. 


QCD Fits Fits to the data of £g. ^ and table show that contributions from second order 
power corrections oc 1/Q^ are not required by the data. Such terms may occur because of 
higher-twist corrections and would be sensitive to the low Q measurements. All event shape 
variables can be well described by first order power corrections oc 1/Q and an exponent p = 0 
for the logarithmic term. 


Note that the power corrections for the jet broadening are expected to have a Q dependence 
different from the other event shapes variables. However, the ansatz oc I/pr fij) 

of eq. (pTf) cannot be supported by the data. The additional logarithmic term varies by a factor 
of ~ 2.7 and leads to a less steep Q dependence, resulting in an unacceptable £t. A modihcation 
of the logarithmic term to ln(/iij//io) with po treated as a free parameter also does not give 
a satisfactory solution. Therefore this logarithmic term is discarded in the following analysis. 
This observation disagrees with the conjectured Q behaviour of ref. 0]. 


The results of the QCD analysis, assuming Pr = Q and /i/ = 2 GeV, are shown in £g. |^. 
The data are well represented by the hts to eqs. (1) - O- The second order perturbative 
QCD predictions are also shown. The power corrections, or hadronisation contributions, are 
substantial at low values of Q, but become less important with increasing energy. A comparison 
with the corresponding LEPTO curves of £g. ^ shows that this event generator adequately 
describes the experimental data. However, the LEPTO parton level distributions look very 
different to the DISENT or MEPJET calculations. 


Theoretical uncertainties in the determination of as{Mz) and do come from the accuracy 
of the QCD calculations and from the choice of different possible scales. Varying the anti¬ 
correlated coefficients Ci and C 2 of eq. within their errors leads to typical changes of = 
±0.001 and = ±0.002. Using a different parton density parametrisation would give 

negligible changes to a* and do. The separation between the non-perturbative and perturbative 
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contributions is characterized by the infrared matching scale ni and the renormalisation scale 
Hr, which should satisfy the relation A /i/ -C Hr- Requiring the criterion hr!Hi > ^ implies 
Hi < 2.5 GeV if hr = Q and hr > 0.8 Q if /i/ = 2 GeV. The upper value of the renormalisation 
scale is chosen such that the htted varies by about the same amount as for the statistical 
error determination, typically hr < 1.5 Q. A variation of the renormalisation scale within 
0.8 Q < Hr < f-f>Q results in shifts of the order of = laool and = ±0.05. The range 
of the infrared matching scale is chosen to be symmetric, Hi = 2.0 ± 0.5 GeV, which leads to 
typical changes of = ±0.002, while the power corrections are proportional to hi- All the 
theoretical uncertainties are added in quadrature. 


5 Results of QCD Analysis 

The hnal results for the power correction parameters do cind the strong coupling constant 
as{Mz) in the MS scheme are compiled in table |^. Note that the overall errors are dominated 
by systematic effects associated with the theory. The largest contribution comes from the 
renormalisation scale dependence, which could be reduced once 0{al) calculations become 
available. However, there remains the interplay between the non-perturbative and perturbative 
regions. An effective means to get a larger separation between the renormalisation and the 
infrared matching scales would be to select data at higher values of Q at the expense of a lower 
sensitivity to both the power correction parameter and the strong coupling constant. 

All event shape analyses give results consistent with each other for do and as{Mz)- It should 
be noted that the htted parameters do are correlated with the preferred value of as{Mz)- The 
correlation coefficients change sign and magnitude for the various event shape variables. 

It is particularly remarkable that the power correction parameters do of the event shapes 
under study are all of the same size to within ±15%. The individual values are consistent with 
being universal and clearly support the concept and computations of ref. 0 of power corrections 
to the mean hadronic event shapes in deep inelastic ep scattering. It is notable, however, that 
the power correction parameter for the jet broadening is slightly lower than those of the other 
event shapes, which may be an indication that the adopted parametrisation is not completely 
correct, as noted in section |. 

The as{Mz) values show a larger spread than expected from the experimental errors alone. 
They become compatible with each other only if the theoretical uncertainties are taken into 
account. This indicates the importance of the influence of higher order QGD corrections, which 
may affect each event shape variable differently. A similar situation was observed in 
event shape analyses performed by experiments at LEP, see e.g. discussion in ref. []^. 

An important result is the 1/Q dependence of the power corrections to the mean event 
shape variables. The similarity of the do values suggests that the assumption of universal power 
corrections may be further tested by performing a common two-parameter £t to both thrust and 
the jet mass data, each of which exploits a different event shape property. The jet broadening 
data are neglected because of the arguments given above. In the common £t no account is taken 
for possible correlations between the event shape variables. The correlation between do and 
as{Mz) is much less than in the individual hts. The common fit (see table |^) gives a relatively 
large y^/ndf ~ 2, a possible reason being the missing higher order QGD corrections. The 
results of the fit are do = 0.491 ±0.003 (exp) IqIoJ® (theory) for the power correction parameter 
and as{Mz) = 0.118 ± 0.001 (exp) ^aooo (theory) for the strong coupling constant in the MS 
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scheme, where the experimental errors are determined from the contour. It is clear from the 
sizes of the experimental and theoretical uncertainties that a better understanding of higher 
order QCD corrections to the event shape variables is necessary for further improvements. 


From hts to the energy dependence {Q = -\/s) of thrust (1 — Tee) and the heavy jet mass 
{Mfj/s) of data 0, and using the the same approach with the QCD prescriptions of 
ref. [ 0 . one obtains values for the power correction parameters which are consistent with the 
ep analyses, see table |^. In a similar analysis with a different choice of e’''e“ experiments 
DELPHI Ends consistent results, but a slightly smaller value of as{Mz) for the thrust 
variable. These observations suggest that the power correction parameters do are universal to 
within ~ 20% in both deep inelastic ep scattering and e+e“ annihilation.]^ 


6 Conclusion 

The event shape variables thrust, jet broadening and jet mass are studied for the hrst time in the 
current hemisphere of the Breit frame in deep inelastic e p scattering. Differential distributions 
and mean values are measured over a wide range of momentum transfers Q from 7 GeV to 
100 GeV. The mean values of the event shapes exhibit a strong Q dependence. They decrease 
with rising Q, i.e. the energy flow in the Breit current hemisphere becomes more collimated. 
The means of thrust and jet mass show a similar energy dependence to that measured in e’''e“ 
annihilation at the scale Q = ^/s^. Differences can be plausibly attributed to the different 
QGD dynamics for, and to the analysis methods which are applied to, ep and scattering. 

Using 0(ag) QGD calculations and a specific model to describe the influence of non- 
perturbative effects on hadronic final states, the Q dependences of the mean event shape pa¬ 
rameters are fitted to determine simultaneously the power or hadronisation corrections and 
the strong coupling constant as{Mz) in the MS scheme without assuming any fragmentation 
model. The power corrections decrease with the expected power of 1/Q and can be described 
by a common universal parameter do which is consistent with the values found in e^e“ experi¬ 
ments. However, the conjectured 1/Q InQ dependence of the jet broadening power corrections 
cannot be supported by the e p data. The present precision of both the power correction para¬ 
meters and the strong coupling constant is limited by theoretical uncertainties due to unknown 
higher order QGD corrections, which are to some extent accounted for by a variation of the 
renormalisation and infrared matching scales. 
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(Q) [GeVI 

1 

1—1 

(Q) |GeV] 

(i-n)/2 

7.47 

0.1678 ±0.0019 ±0.0023 

7.47 

0.2182 ±0.0019 ±0.0024 

8.91 

0.1631 ±0.0021 ±0.0019 

8.91 

0.2009 ±0.0022 ±0.0023 

14.9 

0.1251 ±0.0016 ±0.0021 

14.9 

0.1555 ±0.0018 ±0.0023 

17.8 

0.1193 ±0.0012 ±0.0022 

17.8 

0.1351 ±0.0014 ±0.0025 

24.2 

0.1072 ±0.0012 ±0.0025 

24.2 

0.1125 ±0.0013 ±0.0024 

37.9 

0.0880 ±0.0018 ±0.0021 

37.9 

0.0912 ±0.0019 ±0.0023 

68.0 

0.0746 ±0.0037 ±0.0023 

68.0 

0.0635 ±0.0035 ±0.0049 


(Q) |GeV] 

(Be) 

(Q) [GeVI 

(pc) 

7.47 

0.3566 ±0.0019 ±0.0071 

7.47 

0.1055 ±0.0015 ±0.0023 

8.91 

0.3414 ±0.0023 ±0.0076 

8.91 

0.1009 ±0.0015 ±0.0024 

14.9 

0.2978 ±0.0018 ±0.0069 

14.9 

0.0871 ±0.0010 ±0.0031 

17.8 

0.2704 ±0.0013 ±0.0067 

17.8 

0.0816 ±0.0008 ±0.0020 

24.2 

0.2394 ±0.0014 ±0.0070 

24.2 

0.0780 ±0.0008 ±0.0030 

37.9 

0.2039 ± 0.0023 ± 0.0093 

37.9 

0.0637 ±0.0013 ±0.0030 

68.0 

0.1654 ±0.0046 ±0.0107 

68.0 

0.0512 ±0.0024 ±0.0021 


Table 1: Mean values of the event shape variables thrust (1 — Tc) and (1 — T^)/2, jet broadening 
{Be) and jet mass (pc) as a function of Q. The hrst error is statistical, the second systematic 


Observable 

Cl 

C2 

CLp 

1 

1 — 1 

0.384 ±0.033 

0.57 ±0.21 

1 


0.053 ±0.033 

3.45 ±0.23 

1 

(Be) 

0.990 ±0.121 

2.39 ±0.86 

2 

{pc) 

0.359 ± 0.048 

-0.05 ±0.30 

1/2 


Table 2: Coefficients ci, C 2 and ap of the event shape variables used in the QCD hts 
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Observable 

<yo{pi = 2 GeV) 

ois{Mz) 

xVndf 

HI ep data 

(1 - r.) 

0.497 ±0.005 

0.123 ±0.002 

5.0/5 

{i-n )/2 

0.507 ±0.008 

0.115 ±0.002 

8.5/5 

(Be) 

0.408 ± 0.006 

0.119 ±0.003 

5.3/5 

{pc) 

0.519 ± 0.009 

0.130 ±0.003 

3.1/5 

common £t 
Tc, Pc 

0.491 ± 0.003 

0.118 ±0.001 

39/19 

e+e” data 

(1-Tee) 

0.519 ± 0.009 

0.123 ±0.001 

10.9/14 

{Ml/s) 

0.431 ± 0.020 

0.115 ±0.002 ZZ 

7.8/14 


Table 3: Results on power correction parameters bo and the strong coupling constant as{Mz) 
from fits to the Q dependence of the mean event shape variables. In the common £t possible 
correlations between the event shape variables are not taken into account. The hrst error is 


use the heavy jet mass Mh 


experimental, the second error represents theoretical uncertainties. The e'''e experiments [12 
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Figure 1: Differential event shape distributions of a) thrust 1/N dn/d{l — T^), b) thrust 
1/N dn/d{l — Tz), c) jet broadening 1/N dn/dBc and d) jet mass 1/N dn/dpc- HI DIS ep data 
(full symbols, only statistical errors shown) are compared with LEPTO Monte Carlo simulations 
(—). The spectra for (Q) = 8.3 — 68 GeV are multiplied by factors of 10” (n = 0,4) 
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Figure 2: Mean event shape variables as a function of Q for a) (1 — Tc), b) (1 — Tz)/2, c) 
(Be), and d) (pc)- HI DIS ep data (•, errors include statistics and systematics) are compared 
with LEPTO Monte Carlo simulations for hadrons (—) and partons (■ ■ ■). Data from e’''e“ 
experiments (A) are shown for thrust (1 — Tee)-, calculated for the whole event, and for the 
average of the heavy and light jet masses 
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Figure 3: Differential event shape distributions of a) thrust 1/Ndn/d{l — T^), b) thrust 
1/N dn/d{l — Tz), c) jet broadening 1/N dn/dBc and d) jet mass 1/N dn/dpc- HI DIS ep 
data (full symbols, only statistical errors shown) are compared with second order QCD calcu¬ 
lations (—). The spectra for (Q) = 8.3 — 68 GeV are multiplied by factors of lO” (n = 0,4) 



T> < fl> 






Figure 4: Mean event shape variables as a function of Q for a) (1 — T^), b) (1 — T^)/2, c) 
and d) (pc)- HI DIS ep data (•, errors include statistics and systematics) are compared with 
QCD hts (—) and second order QCD calculations (■ ■ ■) 
















